Crossover in the two-impurity Kondo model induced by direct charge tunneling 
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Quantum critical behavior in the two-impurity Kondo model requires the distinct separation of 
two scales, Tk ^> T* , where Tk is the Kondo temperature and T* is the scale at which the system 
renormalizes away from the quantum critical point to a stable Fermi liquid fixed point. We provide 
a derivation of T* based on the renormalization group to lowest order. This result is confirmed by 
a numerical renormalization group (NRG) analysis which supplements the analytic derivation with 
additional quantitative precision. The form of the low-energy Fermi liquid fixed point is derived 
and subsequently confirmed by the NRG. We discuss implications for series double quantum dot 
systems. 
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I. INTRODUCTION 

Competing orders in correlated electron systems (e.g. 
heavy fermion compounds or high temperature super- 
conductors) lead to new exotic quantum critical points 
(QCPs). Interestingly, QCPs also occur in much simpler 
scenarios where correlations are driven fundamentally by 
a local quantum impurity coupled to a free Fermi sea. 
One of the most fascinating classes of impurity models 
showing quantum criticality are the multichannel Kondo 
models^. At criticality, the electrons are not described 
by a Fermi liquid theory, as has been demonstrated in 
quantum dot experiments 2 . However, one element of 
these critical systems that requires further understanding 
is the nature of the crossover away from the QCP. This 
crossover will generically occur due to arbitrarily small 
symmetry breaking perturbations (such as the presence 
of a magnetic field in the multichannel Kondo effect) at 
temperatures small compared to the crossover scale as- 
sociated with those perturbations. 

Therefore, to obtain closer contact with experiment, 
an accurate estimate of the crossover energy scale is re- 
quired. In this paper, we describe such a calculation for 
a model closely related to multi-channel Kondo models, 
namely, the 2-impurity Kondo model which has possible 
realizations in double quantum dot devices!^ Related 
experiments on quantum dots^ or using a scanning tun- 
neling microscope with one impurity on the tip and one 
on the sample 6 did not observe the physics of the QCP 
so far. Indeed, theoretical predictions of the influence of 
the QCP on the nonlinear conductance^ and shot noise 8 
require the crossover energy to be small relative to the 
Kondo temperature. One of the central purposes of this 
paper is to study the effect of the simplest symmetry 
breaking perturbation in the two impurity Kondo model, 
namely direct charge tunneling, via a numerical calcula- 
tion of the crossover scale. We shall briefly discuss the 
applicability to double quantum dot setups. 



We consider two leads labelled by L and R, each cou- 
pled to one of two spin- 1/2 "impurities" labelled by Sl 
and Sr respectively. The two spins are coupled via a 
Heisenberg exchange interaction K. We also consider a 
direct charge tunneling term Vlr between the two leads. 
The system is described by the model Hamiltonian 



H 



dx ^]^)id x ^/j ni {x) 

j=L,R 

+J(sll-S l + s rr - Sr) (1) 
+KS L ■ S R + Vlr (^,(0)^(0) + ^(0)^(0)) . 

The operator ipjp(x) creates a chiral electron of spin fi 
at position x in the j th lead (these chiral operators are 
defined on the entire real line via the standard "unfold- 
ing" transformation). We have defined the spin operators 

s jj : ~ V^j M (0)— f^jvW (<r are the Pauli matrices) which 
couple to each impurity spin. 

For the special case of Vlr — 0, it is known£ that 
a quantum critical point separating two distinct Fermi 
liquid phases exists for a critical value of K = K c ~ Tk- 
Tk is the Kondo temperature which we take to be 



T K = DV^Je-^r 



(2) 



where 2D is the bandwidth or ultra-violet cutoff and v = 
l/(27r) is the density of states in the leads at the Fermi 
energy (recall that, in our units, the Fermi velocity is 
vp = 1). Throughout this paper, we set h = ks = 
1. For K = K c , the effective low-energy model cannot 
be described as a Fermi liquid. We use the following 
nomenclature to describe the three phases of the model 
at V L r = 0: 

K < K c Kondo screened phase (KSP) 
K = K c quantum critical point (QCP) 
K > K c local singlet phase (LSP) 
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For K close to K c and Vlr close to (the meaning of 
"close" will be defined shortly), the effective description 
of the system will depend on the temperature T. For 
T Tk, the spins are only weakly coupled to the leads. 
For f«r« Tk, the system will be described by the 
QCP while, for T <C T* , the system will be described by 
a Fermi liquid that is continuously related to the KSP or 
the LSP via the parameter Vlr as will be described in 
detail in § HIl 

It has been proposed^ that the temperature scale T* 
can be estimated as 

T* = a {K ~ Kc)2 + bT K {vVLR? (3) 

where we have inserted two dimensionless numbers a and 
b that are expected to be of order unity. The value of 
T* determines how close K must be to K c and how close 
Vlr must be to in order to observe the QCP: the values 
must be such that T* < T K . 

To understand where this estimate for T* comes from, 
let us assume that K is tuned close to K c and that vVlr 
is very small. We want to estimate the crossover en- 
ergy scale at which the system renormalizes away from 
the critical point, due to both SK = K — K c and also 
vVlr- To make a rough estimate, we use the lowest or- 
der weak coupling renormalization group (RG) at energy 
scales above Tk and the RG at the QCP below Tk- At 
weak coupling, 5K has dimension 1 and vVlr is dimen- 
sionless. Thus, at scale Tk, 

SK(T K ) « (D/T K )SK 
vV L r{T k ) « Wlr. (4) 

Here, the quantities on the right are the bare ones at scale 
D and the quantities on the left are the renormalized 
ones at scale Tk- At energy scales E < Tk, both 5K 
and vVlr have dimension 1/2 so 

SK(E) « (D/TkXTk/E^SK 
vV L r{E) « {Tk/EY^vVlr. (5) 

We estimate the crossover energy scale T§k associated 
with the coupling SK (not to be confused with the Kondo 
temperature Tk) by setting 5K(Tsk)/D « 1. Similarly, 
we estimate the crossover energy scale Tlr associated 
with the potential scattering by setting vVlr{Tlr) ~ 1. 
Thus 

T SK « a{5K) 2 /T K 

T L r « KvVlr) 2 T k (6) 

where the dimensionless factors a and b are used to 
account for quantitative details missed by this simpli- 
fied analysis. The total crossover scale is then given by 
T* = Tsk + T L r. We might expect further modifica- 
tions if we took into account higher order terms in the 
RG equations at weak coupling and in the vicinity of the 
QCP. Hopefully this just leads to corrections which are 
logarithmic in the dimensionless parameters. 



The estimate of eq. Q agrees with similar analytic 
arguments^iS (although it is slightly smaller than that 
found in the latter reference). In §|H]we verify this esti- 
mate by a transparent and systematic NRG study with 
the aim of estimating this crossover scale. We find sat- 
isfactory agreement with eq. (j3]) though with the mag- 
nitude of the coefficients a and b differing from unity. 
We also find evidence that the coefficient b may have 
some residual dependence on the Kondo coupling J that 
is not explained by the simple scaling analysis that lead 
to eq. ©. 

In § IIII1 we derive the Fermi liquid theory for the stable 
low-energy fixed point which is parametrized by the value 
of Vlr- The resulting theory is described in terms of a 
boundary condition at the origin which can be related to 
a scattering phase shift. The formula for this phase shift 
is predicted analytically and supported by comparison 
with the phase shift derived from the NRG. In § IIVI wc 
discuss the two impurity Anderson model version of this 
problem, showing that the simplified Kondo model of Eq. 
@ is not the correct low energy approximation to the 
Anderson model. 



II. ESTIMATING T* USING THE NRG 

The numerical renormalization group is a powerful, 
non-perturbative algorithm for studying quantum impu- 
rity systems. It was originally developed to study single- 
impurity model in [13l4l5l | where the technique is exhaus- 
tively described. It has since been applied to numerous 
other impurity problems 1 - including the two-impurity 
problem that we study tiered. We refer the reader to 
these references for details on the implementation of the 
NRG and only review those elements necessary for cal- 
culating T*. 

The key idea of the NRG is to approximate the original 
Hamiltonian describing the leads by that of two semi- 
infinite tight-binding 'Wilson chains" 

H 1 °° 

^ = 2( 1+A_1 )E E A-f£n(/4j»+i )W +b.-c..) 

n— Op— e,o 

+ vJ E /oV^r/o P „ • (S L + S R ) 

p— e,o 

+"J (fL^foou + h.c.) • (S L - Sr) 

+ ~D Sl ' Sr + 2vVlr (aW-W - fLfifoof^ ■ (7) 

Each of the / npM operators are complicated linear com- 
binations of i/v ~ tyitu + V'au for p = e and -0 OM ~ 
ipLfi — ipRfj, for p — o (see ref. |9| for more details). Here, 
2D is the bandwidth cutoff, A > 1 is a numerical pa- 
rameter (we use A = 3), and ~ 1 is a dimensionless 
function of A. The limit A — » 1 recovers the original 
model, albeit in a discrete basis. 

To simulate the RG flow numerically, one introduces 
a series of dimensionless Hamiltonians by truncating the 
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semi-infinite chains to TV + 1 sites each and setting the 
overall energy scale to be of order unity: 



'jV-1 



Hn = A^)/ 2 £ £ A"**, (fl m U + i, m + h.c.) 

^ n— p=e,o 

p— e,o 

+ ^ (/ol^/oo, + h.C.) • (S L - S R ) 

+KS L ■ S R + V LR (fljoe^ - /oViW) } ■ (8) 

The parameters J, K, and Vlr are rescaled versions of 
J, K , and Vl_r such that 



— = lim — 

D Af^oo 2 



l + A^lA-^- 1 )/ 2 ^ 



N- 



(9) 



The RG transformation is realized by numerically di- 
agonalizing the series of Hamiltonians Hn, starting with 
N = 0, using the lowest^ eigenvalues (and associated 
eigenvectors) of Hn to define the Hn+\ matrix, until 
the energy spectrum ceases to change from one iteration 
to the nexfei£. The energy scale of the spectrum of the 
dimensionless Hn is of order unity, meaning that it de- 
scribes H at a temperature scale Tn given by eq. ([9]) to 
be 

T N ~ ^(l + A-^A-^- 1 )/ 2 ^. (10) 

This identification will allow us to use the flow of the 
NRG energy levels to measure the energy scale T*. 
Herein, we will take D — 1 and measure all energy quan- 
tities in units of D. 

To see how this works in practice, consider the flow 
of energy levels when the Hamiltonian parameters have 
been tuned close to the quantum critical point. A few 
such energy levels have been plotted in Fig. [TJ The en- 
ergy levels flow close to those of the unstable QCP but 
eventually flow to a stable fixed point (for Vlr = it is 
either the LSP or KSP depending on the initial value of 
K; see black solid lines) at an iteration number TV*. We 
identify the KSP and the LSP by looking at the quan- 
tum numbers and degeneracy of the energy levels in each 
fixed point regime and comparing to those expected from 
the Wilson chain forms of the fixed point Hamiltonians 9 . 
This is not possible for the QCP which does not allow 
a simple Hamiltonian form. Identifying the spectrum of 
the QCP requires more sophisticated techniques^ 9 -. 
Quantitatively, we measure the value of TV* when the 
energy level crosses over from its QCP value to its 
stable fixed point value, then take the average (TV*) of 
the values TV* for the first 20 energy levels. The resultant 
value of T* is determined from eq. (ITU|) to be 



^(l + A-^A-^- 1 )/ 2 



(11) 



To estimate the value of a in eq. ((31), we set Vlr — so 
that T* = Tsk and choose a value of v J. We then obtain 




-r| 2c (0) 




FIG. 1. The lowest three NRG energy levels in the charge— 1-1, 
spin- 1/2, parity-even subspace as a function of even iteration 
parameter TV". A value of vVlr = was used to determine 
the solid black lines whereas a value of vVlr = 3.372 x 1CP 7 
was used for the red dotted lines. The arrows on the left axis 
label the prediction of the second and third excitation of the 
QCP from ref. Qjl (the first excitation is used to match the 
overall scale of the predicted spectrum with that from the 
NRG). The arrows on the right axis label the first and second 
single-particle excitations used in the determination of the 
even channel phase shift, eq. (|21fl . The arrows on the bottom 
axis label the value of (TV*) used to estimate T* for each of the 
two spectra. Here, vj = 0.15 and K = 1.30096478 x 10" 3 > 
K c so that the final, stable fixed point is that of the LSP. 



NRG spectra for a series of Hamiltonians with differing 
values of K, starting with 20 K = K c then tuning K 
away from K c until T* ~ Tk [i.e. when the QCP is 
no longer reached in the RG). From these spectra, T* 
is extracted as described above and a plot is made of 

InT* vs. In (K - K c f /T K ; see figure H A linear fit 

is made to the data and compared with the expectation 
from eq. jSJ): 



InT* 



m„ In 



{K-K c 



In a, V] 



LR 







(12) 



where Tk is determined from the input value of v J us- 
ing eq. ((2|). Obtaining the slope m a — 1 from the fit 
provides a check on the dependence of T* on K — K c 
while the intercept of the fit provides an estimate of the 
dimensionless coefficient a. 

A similar procedure is used to measure b. With 
K = K c so that T* = T LR , NRG data is obtained for 
a series of values of Vlr, starting with Vlr = and in- 
creasing Vlr until T* ~ T K - A plot is made of InT* 
vs. In [Tk(vVlr) 2 ] ; see figure 03 Again, a linear fit is 
made to the data and compared with the expectation 
from eq. ([3]): 

InT* = m 6 ln [T K {vV lr ) 2 ] + ln6, K = K c . (13) 

As before, obtaining the slope mi, — 1 provides a check 
on the dependence of T* on vVlr and the intercept de- 
termines the value of b. Consistent values of a and b over 
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ln[(K-K)VT K ] '"[ T k( vV lr>~J 



FIG. 2. Data for T* as determined from the NRG with Vlr = 
0. The solid line is the best linear fit to all of the data. 




ln[T K (W LR )-] 



FIG. 3. Data for T* as determined from the NRG with K = 
K c . The solid line is the best linear fit to all of the data. 



several values of v J provide confirmation of the proposed 
crossover energy scale Q . 

We have carried out the procedure described above for 
four values of v J. The parameters a and b are tabulated 
in Table U together with the error arising from the fit to 
the data. In Figures [5] and [3] we have plotted the data 
for all four of these iterations and find reasonably good 
data collapse. 

Before describing our analysis of T* , it is interesting 
to note that, although we find that K c ~ Tk as is often 
quoted in the literature, we do not find the consistent 
value of K c = 2.2Tr- as reported in ref. [jj Rather, we 
find K c — a(vJ)Tji where a ranges monotonically (at 
least according to the four values obtained and listed in 
Table U from a(0.217) = 3.36 to a(0.100) = 1.8. Al- 
though this does not change our conclusions regarding 
T* , we simply point out that the relation between K c 



FIG. 4. An enlarged version of a section of Figure [3] 

and Tk may not be as simple as first presented^. 

Returning to our analysis of T*, it is seen that, while 
the value of a is roughly of order unity, the value of b 
here is two orders of magnitude larger than unity. Fur- 
thermore, while the values of a appear to be consistent 
for all values of v J, there seems to be a slight trend of b 
decreasing with decreasing J. This can be seen in Fig- 
ure [4] which is simply an enlarged area of Figure [3] The 
fact that we consistently obtain nib ~ 1 means that the 
T* dependence on vVlr is certainly quadratic but that 
the J dependence of the coefficient in the second term of 
eq. ([3]) may be different than that included in the factor 
of Tk- 

For further analysis, we look at the explicit J depen- 
dence of the coefficient of {vVlr) 2 in eq. (J3]) to see how 
well it matches that of the predicted VTk- To do this, we 
set K = K c and plot in Figure[5]the value of T* / (vVlr) 2 
extrapolated to vVlr = for each of the four data sets 
versus the corresponding value of To this data we 

have fit a function of the form 

2 = bT K = bV^Je-^r (14) 

{vVlr) 

and obtained a value of b — 118. This value is close to the 
mean of the four values of b listed in TableUand the above 
function provides a reasonable fit to the data. However, 
as before, there is a hint of further J dependence as the 
function appears to overestimate the vJ = 0.1 data while 
underestimating the data at vJ — 0.217. 

III. FERMI LIQUID THEORY 

Having determined the two energy scales Tk and T* , 
the latter with aid from the NRG, one can perform per- 
turbative calculations in wide temperature ranges using 
effective theories appropriate for each regime. In the 
high temperature weak coupling regime T 3> Tk, one 
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vJ 


T K (xl0~ 3 ) 


K c (xl0~ 3 ) 


m a 


Am 


a 


Aa 


mj, 


Arab 


b 


Ab 


0.217 


4.644 


15.6194231 


1.004 


0.004 


0.35 


0.04 


0.994 


0.001 


130 


6 


0.183 


1.811 


5.55847415 


0.999 


0.005 


0.33 


0.05 


0.997 


0.002 


125 


8 


0.150 


0.4929 


1.30096469 


1.001 


0.004 


0.37 


0.04 


0.999 


0.001 


114 


3 


0.100 


0.01436 


0.02617379 


0.99 


0.01 


0.4 


0.2 


1.009 


0.004 


100 
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TABLE I. NRG results for four values of uj. Note that all temperatures are measured in units of D with 2D being the 
bandwidth in the leads. The dimensionless parameters a and b are defined in eq. (0 and Aa: indicates the error in x arising 
from the linear regression. 




FIG. 5. The value of T*/(vV L r) 2 extrapolated to uVlr = 
for each of the four data sets (points). Here, K = K c . A 
function of the form of eq. (|14|l is fit to the data (solid line) 
with the free parameter found to be b = 118. 



can apply perturbation theory in J to the original model 
described by the Hamiltonian of eq. (1). In the interme- 
diate temperature QCP, T*«T« Tk, one should use 
the conformal field theory^ 2 - (or the corresponding the- 
ory in the language of abelian bosonization 23 ) and apply 
perturbation theory in the irrelevant operator. In the 
low temperature Fermi liquid (FL) phase, T < T*, one 
can use a Fermi liquid theory in terms of single fermion 
scattering states as briefly reported in ref. 

In this section we shall give a derivation of the ef- 
fective FL Hamiltonian (see eq. (f5Tj)). In sub-section 
A we derive the phase shift which characterizes the 
Fermi liquid fixed point, and which depends on the ratio 
(K — K c ) I \TkvVlr\ ■ This is done using only the formula 
for the T = conductance on the line of Fermi liquid 
fixed points derived in ref. Q, a standard Fermi liquid 
formula for the conductance in terms of the phase shift, 
and a symmetry argument. In sub-section B, we derive 
the leading irrelevant interactions and corresponding cou- 
pling constants on the line of Fermi liquid fixed points. 
This is done starting with the abelian bosonization de- 
scription of the QCP— ^ which uses Majorana fermions 
that are non-locally related to the original Dirac fermion 
fields in eq. ©. We will then see that, at temperatures 



lower than T* and with the relevant operator present, 
single electron scattering states become the correct par- 
ticles in terms of which the FL theory is conveniently 
written. 



A. Phase Shift Analysis 

Before describing the derivation of the FL Hamilto- 
nian, we present a simple derivation of the form of the 
phase shifts in terms of which the single electron scatter- 
ing states are defined. These are then compared to the 
phase shifts extracted from the fixed point NRG spec- 
trum. Since one can derive the same form for these phase 
shifts from the fixed point analysis described at the end of 
this section, numerical confirmation of the phases shifts 
provides additional support for the analytic calculation 
of the FL Hamiltonian. 

In the simpler FL theory of the single channel Kondo 
effect^, the Hamiltonian is written in terms of weakly 
interacting fermionic scattering states which are simply 
the original Dirac electrons in which the zero temper- 
ature scattering phase shift is incorporated. The same 
holds true for the two-impurity model under consider- 
ation. Using the L R symmetry, the original Dirac 
fermions satisfy a FL boundary condition (BC) 

(^ + ^r m )(0 + ) = e 2iS '(^ Lll + V^)(0-), 

(V^ " ^m)(0 + ) = e M °(^L M - ^)(0")- (15) 

Furthermore, using the special particle-hole symmetry 

"> V'RA. -> (16) 

it follows^ that S e = — 5 = 5. 

To calculate the phase shift 8, we use its relation to 
the zero temperature conductance^ 

A G = sin 2 («s e _ So ) = sin 2 25. (17) 
Comparing this withi 
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we can immediately extract the form of 5 
S = 2 ai 'g<V T <5if + W t lr) 



-arg 




wVlr 



(19) 



where, in the last equality, we have substituted the ex- 
pressions for Tgx and Tlr, from eq. 

For Vlr > 0, the phase shift 5 changes from to ir/2 as 
function of K, and it takes the value of ir/4 at K = K c . 
This agrees with the numerical results of Jones el. alM 
While the original electrons suffer a phase shift, the single 
particle scattering states incoming from lead j = 
1,2 = (L,R), defined by 



cos 25 
—i sin 25 



-i sin 2(5 
cos 25 , 



are continuous at the origin: ^>i ll (0 + ) — ^^(0 ). 



(20) 



The form of this predicted phase shift, eq. (|T9|i . can 
be compared with the phase shift derived from the NRG 
fixed point spectrum. To extract the phase shift, one 
looks at the many-body NRG energy spectrum for the 
Hamiltonian Hn of eq. ((8]) where N is large enough such 
that the RG has reached one of the FL fixed points. Un- 
like the QCP, the LSP and KSP many-body spectra are 
made up of two channels of single-particle/hole excita- 
tions combined in such a way so as to respect Fermi 
statistics. By looking at the quantum numbers of the 
lowest many-body NRG energy levels, one can determine 
the two lowest single-particle/hole excitations r\\ v and i]2 P 
in each of the p = e (even), o (odd) channels. From these 
we define the phase shift as 



5 P = Q 



r)i P + r)2 P 



(21) 



for the case of even N. Here, Qi p = ±1 is the charge of 
the lowest excitation in the p channel indicating whether 
the spectrum is shifted up (Q\ p = 1, 6 P > 0) or down 
(Qip = — 1, S p < 0) relative to the spectrum with Vlr = 
0. The ry's are marked in figure [TJ The definition for odd 
N is simply shifted by ir/2. This definition of the phase 
shift at the fixed point follows closely that used in ref. [2?] 
for a model of a quantum dot in an Aharonov-Bohm ring. 

In Figure [HI we have plotted the phase shifts as derived 
in this way from the NRG and compared them with those 
predicted from eq. (TT91 for both K > K c and K < K c . 
We find the agreement to be quite good and take this as 
support for our analysis of the FL fixed point. 

Slight deviations from the continuum model used to 
derive eq. (TT91 are known to exist due to the fact that 
A > 1 in the discrete NRG Hamiltonian (see ref. [27] for a 
discussion of this effect). Furthermore, we obtain better 
agreement with the NRG for the phase shifts in the LSP 
than we do for phase shifts in the KSP. The nature of this 
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FIG. 6. The phase shifts S e and S in the even and odd 
channels (respectively) of the FL fixed point as derived from 
the NRG (points) with vj — 0.15 as well as those phase 
shifts predicted analytically by eqs. (|15|) and ()19[) (lines). The 
values for Tk, K c , a, and b appearing in eq. (|19p are taken 
from the v. J = 0.15 line of Tabled] We have done this for 
both the LSP, K = 1.30096478 x 10~ 3 > K c , and the KSP, 
K = 1.30096452 x 10~ 3 < K c ). 



discrepancy looks very similar to that due to the presence 
of an additional potential scattering term that is gener- 
ated by the Kondo interaction in the screening channel 
when particle-hole symmetry is broke n 27 ' 28 . Since the 
presence of a non-zero Vlr breaks particle hole symme- 
try, one would expect such an effect but only in the KSP 
where Kondo screening occurs. This is precisely what is 
seen in Figure [HI However, the form of this additional 
potential scattering was only derived for single-channel, 
single-impurity models 27,28 so more analysis is required 
to determine for certain if this is the nature of the dis- 
crepancy in the KSP phase shifts. Nevertheless, it is clear 
that eq. (fl"9]) captures the leading order contribution to 
the phase shift for the entire manifold of fixed points. 



B. Fermi Liquid Hamiltonian 

We now turn our attention to the derivation of the 
FL Hamiltonian. Using abelian bosonization, one can 
write the original free Fermion theory (eq. ([2|) with 
J —> and Vlr — > 0) in terms of 8 chiral Ma- 
jorana Fermions xt associated with the real (xi = 

) an< ^ imaginary (xt — ~ ) P ar t s of the 

charge, spin, flavor and spin-flavor fermions ipA ~ 
e~ z ^ A , (A = c,s,f,X). The bosonic fields cf>A are lin- 
ear combinations of the four bosonic fields associated 
with the original Dirac fermions, ipip ~ e -i( ^, (i = 
L,R= 1,2, fx =t,4= 1,2) given by {0 C > S ,0/,0 X } = 

The free Hamiltonian is iJo[{x'}] = 
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{xf ,Xi,X2>xf )Xi»X2>Xi)X2} is an arbitrary rela- 
beling of the 8 fields. Turning on the Kondo coupling 
J at Vlr = 0, the QCP is obtained at K = K c . It is 
described simply in terms of a change in the BC relative 
to the free case in which Xi(0~) = Xi(0 + )i — 1; — 8)- 
The change in BC occurs only for the first Majorana 
fermion, x'i(0 _ ) = — xi(0 + ). In our chirality convention 
for the 1 dimensional fields, the region x > (x < 0) 
corresponds to the incoming (outgoing) part of the field. 
For energies <C Tk, we define a new basis 

Xi(x) = x'i(x) sgn(x), 

Xi(x) = Xi(x), (z = 2, . . . , 8) (22) 

and write the leading terms in the Hamiltonian describ- 
ing deviations from the QCP due to finite K — K c as well 
as finite Vlr as H QC p = H [{x}] + SH QC p with 7 

2 

6H QCP =i z 2 X iXi(0)a. (23) 

i=l 

Here, a is a local Majorana fermion, a 2 = 1/2. Both 
terms in SHqcp have critical dimension 1/2 so that they 
destabilize the QCP. The coupling constants satisfy Ai oc 
K — K Cl \ 2 oc Vlr and 

T SK = Xf, T LR = Xl T* = X\ + \\ (24) 

where T* is also given in cq. ([3]) with the coefficients 
a and b determined numerically in Table 1. Below the 
crossover scale T* , the system flows to FL fixed points 
whose nature depend on the ratio A1/A2. 

The crucial observation is that only the linear combina- 
tion [XiXi(x) + X 2 x 2 (x)}/ X of the 8 Majorana fermions at 
the quantum critical p oint participates in this crossover. 
Here A = yfif+Jg. It can be shown that the ef- 
fect of SHqcp is to modify the BC for only this linear 
combination by a simple sign change at the boundary, 
(Aixi(0+) + A 2 x 2 (0+))/A = - (AiXi(O-) + A a x 2 (0-))/A. 
In order to write down the FL fixed point Hamiltonian, 
we define a new basis with modified BC, {rj}, where 

r)i(x) = sgn(a;)(AiXi(x) + A 2 X2(a;))/A , 

= {-X 2 xi{x) + XxX2{x))/X, 
T H (x)= Xi (x), (i = 3,...,8). (25) 

We can therefore write the Hamiltonian for the FL fixed 
points as Hfl = Ho [{f}}] + SHfl- For clarity, we re- 
capitulate the notation: in the weak coupling regime we 
defined the set of Majorana fermions {x'}', in the QCP, 
the Hamiltonian eq. (|23| is written in terms of {%} given 
by the relation eq. (f2"2"| ; in the FL regime, the Hamilto- 
nian will be written in terms of the {77} operators defined 
in eq. (J2SJ). 

Near a FL fixed point the leading interactions have 
scaling dimension 2, which lead to corrections of various 
quantities which vanish at zero temperature as T . In 
our local theory, the interaction SHpl acts only at x = 



and involves uniquely rji which is the only field partic- 
ipating in the crossover in eq. (|2"3"]) . The only possible 
such operator is 

SH FL = X FL imdxm \x=o- (26) 

In the FL theory, T* acts as a high energy cutoff. From 
dimensional analysis, the coupling constant scales as 
Xpl ~ w- Since eq. (|23p is quadratic, we can deter- 
mine the coefficient exactly by matching at low temper- 
atures the results of calculations of a physical quantity 
(e.g. the conductance of a quantum dot 7 -) as calculated 
either using eq. (|2"3"|) or eq. ([2T5)) . In this way, we obtain 
Xpl = y^- 

We have derived the FL scattering states ^j^, j = 
1,2 = L, R, in eq. pt)]) and now rewrite the FL interac- 
tion of eq. (|26l) in terms of them. First we write SHpl 
using eqs. (f22|) and ([25]) as 

SH FL = ^[cos 2 (2S)(x^d x x / 1 )+sm 2 (2S)( X ' 2 td x x , 2 ) 

+ sgn(x)^ sin(45)(xiidxX2 + X^xX'i)]- (27) 

These quadratic operators of dimension 4 can be written 
as a product of two normal ordered quadratic forms using 

XiiOxXi + X 2 iOxX 2 = (« : XiX 2 > 

X'iid x x'i ~ X^dxX'2 = (i ■ Xixs O 2 ~ (i : X2X3 0", 

Xi^xX 2 + X 2 ^xX'i = {: X1X3 =, = X3X2 ■}+, (28) 

where {, } + stands for the anticommutator. 

These normal ordered quadratic forms are related to 
the flavor currents: 

■h X : Kf*,> ■■ (29) 

where Pauli matrices. Indeed in the range x > 

at which the scattering states coincide with the original 
fermions *S?(x > 0) = ip(x > 0), these can be expressed 
as 

J]=i- X1X2 J} = i: XW 3 Jf = i • X3X2 : ■ (30) 
Using eqs. ([27). (l28t and ([5D]). we finally obtain 

SH FL = 7^ffMJ f , (31) 

where 

(l \ 

M = cos4<5 -sin 45 • 
y0 -sin 4<5 -cos4<5y 

One can see that, due to the potential scattering term 
oc VlrJj, the flavor SU(2) symmetry is reduced down 
to U(l). One can rewrite SHpl in different ways, for 
example, in ref. 0, the same equation is written in an 
explictly spin SU(2) symmetric way. 
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We emphasize the universality of the derived FL 
Hamiltonian: for any value of the ratio of original pa- 
rameters Vlr/{K — K c ), all coupling constants of Hrl 
are determined up to the overall energy scale T* which 
was calculated numerically here. The universality fol- 
lows from strong restrictions due to a large symmetry 
that emerges close to the quantum critical point^ 2 - and 
leads to the simple form of Hqcp in eq. (1231) . One can 
see that Hq p has an SO (8) symmetry. Due to SHqcp, 
the crossover from the QCP to FL fixed points has an 
SO(7) symmetry represented by rotations of the vector 
(X2)---)Xs)- This symmetry considerably restricts the 
possible interactions and sets relations between the dif- 
ferent coefficients in eq. (|3"Tj) . The conditions of validy of 
5Hfl is small deviations from the QCP \K — K c \ <C Tr- , 
vVlr *C 1. In addition, a scale separation T* <C Tk is 
required. This scale separation was shown to hold using 
our NRG calculations as discussed in detail in ^ HTl 

In practice, the SO (8) symmetry at the quantum crit- 
ical point is broken by marginal and irrelevant opera- 
tors at the QCP such as the leading irrelevant operator 
(dxXi) a at the quantum critical point. However, these 
will be associated with a small parameter 1/Tk and, 
hence, are neglected at T* <C Tk- For finite Vlr or 
K —K c , additional marginal and irrelevant terms are pro- 
duced at the quantum critical point, part of which were 
present before. However, close enough to the quantum 
critical point, namely for vVlr <C 1 and \K — K c \ <C Tk, 
those perturbations can be safely ignored. 



IV. THE TWO IMPURITY ANDERSON MODEL 

We have presented a detailed analytic and numerical 
calculation of the scale for the crossover away from the 
QCP in the two-impurity Kondo model. It is at this scale 
that the system flows to a Fermi liquid fixed point on 
a one-dimensional manifold of fixed points parametrized 
by Vlr/ (K — K c ). The Fermi liquid behavior along this 
manifold has been derived in detail and confirmed nu- 
merically. Furthermore, the more precise calculation of 
T* presented in this paper provides a more accurate esti- 
mate of the linear conductance written explicitly in terms 
of the crossover scales T* and Tlr, which in the limit 
T < T* is given byi 



G 



2e 2 T LR 
h T* 



( 2nT 



(32) 



As discussed in the previous section, the T = limit of 
this expression for the conductance agrees well with our 
NRG results. 

We now discuss the applicability of our analysis to 
double quantum dot setups. First, real quantum dots 
include charge fluctuations and hence are described by 
the Anderson model rather than Kondo model. Second, 
in real quantum dots there are additional tunneling pro- 
cess not included in our model with the direct tunnel- 



ing Vlr term. In those additional terms the tunneling 
electron's spin interacts with the impurities. We will re- 
express our Kondo model results in terms of Anderson 
model parameters^ We define the Anderson model with 
hopping ta between the conduction electrons and the im- 
purities and hopping tLR between the two impurity sites 
on which there is a Coulomb repulsion U and resonant 
energy level Ed = —17/2. Carrying out a Schrieffer- Wolff 
calculation up to third order in the tunneling terms, the 
Kondo parameters are given by: 



J 
K 



8t 2 d /U 

<r/U 





Tj 



K 



U 



■KVj 



exp 



1 

uJ 



(33) 



(See ref. [23 for the final formula in Eq. (f3"3"|) which is the 
same as that used in ref. [l2|.) Interestingly, the direct 
tunneling operator is completely absent at Ed = — U/2. 
Instead, one obtains the tunneling terms 



5H = V&Wl^ + ^lJSl 



Si 



(34) 



+ VLRi-^L^^Rp-' +i^R^^L^,') ■ (Sr xSl), 
describing tunneling mediated by interaction with the im- 
purity spins. Here V? R = V^r = 1G % LR ■ 

We now analyze these operators from a field theoret- 
ical perspective, assuming that originally the coupling 
constants are very small, and check if they eventually 
destabilize the fixed point by turning on the relevant di- 
mension 1/2 operator with coupling constant A2. It is 
convenient to introduce three types of particle hole sym- 
metry: 



phi 
ph2 
ph3 



l/) -> (T V T Z lp\ 



(35) 



where the first two were introduced in Ref. [T^ [here 3(f) 
are Pauli matrices acting in spin (channel) space] . Both 
terms in SH are odd under phi and even under ph2. 
The relevant tunneling operator with coefficient A2, has 
the same symmetry as Vlr^ which is also odd and even 
with respect to phi and ph2, respectively. To distinguish 
between the two terms in 5H we introduce ph3. The 
operators with coupling constants Vlr and V^ R are both 
odd under p-h3. The relevant tunneling operator with 
coefficient A2, has the same symmetry as Vlr^ and must 
therefore also be odd under p-h3. Hence V^r should 
generate A2 . It can be checked that the V^ R operator is 
even under p-h3; hence it can not generate A2. Since both 
coefficients V^r and V^ 1 are initially small, we keep 
only the former which will grow under renormalization 
group near the critical point. Naively, we might expect 
that we could replace Sl ■ Sr by its expectation value at 
the quantum critical point in calculating A2: Sl • Sr ~ 
— 1/4, corresponding to equal probability of the spins to 
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be found in a singlet or in a triplet stated Hence, SH 
effectively becomes 

SH -> + (36) 

Based on this naive order of magnitude estimate it is in- 
teresting to discuss how our results for T* affect the feasi- 
bility of observing the QCP in such a quantum dot exper- 
iment. The criterion for the observability of the QCP is 
the separation of the two crossovers, namely the crossover 
from weak coupling to QCP occuring at scale Tk and the 
crossover from QCP to FL occuring at scale T* . This re- 
quirement reads Tk/T* 3> 1. Using the Kondo- Anderson 
correspondence with Eq. (|36|) . the crossover temperature 
is: 

_ bT K {vJf 

R ~ — up — LR ' ^ ' 

Evaluated at ti,R — t c LR = yJUK c /A and estimating 
K c = 4(tl R ) 2 /U pa 3T K and b pa 1 15, we obtain 

^ 20 ZM^)I. (38 ) 



In order to estimate the magnitude Tk/T*, we use a 
typical value of U = 1.5meV and a value of vJ = 0.217. 
Rather than estimating Tk using the Kondo model for- 
mula dU , we use the Kondo temperature of eq. (|33|) as de- 
rived from an Anderson model 29 which is more applicable 
to quantum dots. In this case, we obtain Tk « 3.1/jeV pa 
0.04K so that the ratio of scales is T K /T* pa 500 » 1. Al- 
though it would be very difficult to attain temperatures 
T < Tk for such a small value of the Kondo temperature, 
holding U — 1.5meV fixed in eq. {55} suggests that one 
can obtain a ratio of Tk/T* sa 10 at a Kondo tempera- 
ture of T K ~ iO^eV pa 0.5K (VJ ?a 0.4), a temperature 



that can be obtained in modern quantum dot experi- 
ments. Hence, we would conclude that it is possible to 
realize the QCP in a double quantum dot system. 

However, NRG calculation !] 11 ' 12 ' 31 on the two impu- 
rity Anderson model show that in order to achieve scale 
separation Tk ^> T* , much smaller values of vJ oc Y/U 
are needed in the Anderson model; whereas in the for- 
mer paragraph it was estimated that vJ — 0.4 would be 
sufficient to fulfill a scale separation, it appears that vJ 
should be at least one order of magnitude smaller, cor- 
responding to Y/U — 0.05 in the Anderson model^i to 
achieve such scale separation. 

This means that the operator SH in Eq. (f34|) turns 
out to lead to a much larger crossover energy scale T* 
as compared to the direct tunneling Vlr term. Thus, we 
conclude that in order to observe the quantum critical 
point in quantum dots one should achieve the weak tun- 
neling and very low temperature regime, which goes be- 
yond the current ability. Otherwise, another parameter 
must be fine-tuned in order to make A2 small. Alterna- 
tively, more elaborate schemes may be used, where the 
two leads are separated by multiple quantum dots hence 
suppressing the tunneling. 4 Further analysis is required 
to explain why the operator SH in Eq. (|34[) above and 
the Vlr term in Eq. (1) give so different crossover energy 
scales. 
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